Effects of Resonant Cavity on Macroscopic Quantum Tunneling 
of Fluxon in Long Josephson Junctions 
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We investigate the effects of high-Q c resonant cavity on macroscopic quantum tunneling (MQT) 
of fluxon both from a metastable state to continuum and from one degenerate ground-state of a 
double-well potential to the other. By using a set of two coupled perturbed sine-Gordon equations, 
we describe the tunneling processes in linear long Josephson junctions (LJJs) and find that MQT in 
the resonant cavity increases due to potential renomalization, induced by the interaction between 
the fluxon and cavity. Enhancement of the MQT rate in the weak-coupling regime is estimated by 
using the experimantally accessible range of the model parameters. The tunneling rate from the 
metastable state is found to increase weakly with increasing junction-cavity interaction strength. 
However, the energy splitting between the two degenerate ground-states of the double- well potential 
increases significantly with increasing both the interaction strength and frequency of the resonant 
cavity mode. Finally, we discuss how the resonant cavity may be used to tune the property of 
Josephson vortex quantum bits. 

PACS numbers: 74.50.+r, 74.78.Na, 85.25.Cp 



I. INTRODUCTION 

Experimentally observed 2 quantum behavior of 
Josephson vortices (i.e., fluxons) at ultra-low tempera- 
tures has opened up a possibility of realizing quantum 
computers based on long Josephson junctions (LJJs). 
This observation led to much interest on Josephson 
vortex quantum bit* 3 - - — (qubit) as an alternative to the 
previously proposed superconducting qubits. Similar 
to other approaches based on Josephson junctions such 
as charge^ phased and flux 8 - qubits, Josephson vortex 
qubit (JVQ) is also a promising candidate for quantum 
computation application. Due to its weak interaction 
with decoherence sources in the environment at low 
temperatures, the JVQ may have significant advantages 
over the other superconducting qubits. For instance, a 
significantly longer decoherence time was suggested as 
one such advantaged 

The JVQ takes advantage of the coherent superposi- 
tion of two spatially separated states arising from the low 
temperture property of a trapped fluxon in a double- well 
potential. This property includes (i) energy quantiza- 
tion and (ii) macroscopic quantum tunneling 2 (MQT) . We 
note that, for linear LJJs, the fluxon potential for either 
metastable state or JVQ may be obtained 9 - by using Nb- 
AlO^-Nb junctions and by implanting either one or two 
microresistors in the insulator layer, respectively. For ap- 
plication of JVQs, tuning both the decoherence time and 
the level of entaglement by controlling the qubit prop- 
erty is essential. However, due to its weak interaction 
with external perturbations, an effective tuning mecha- 
nism for JVQ is less clear. Recent studie s 10 ' 11 on using 
microwave cavity for both tuning a single phase qubit 
and inducing interaction between either two charge or 
two phase qubits suggest that resonant cavity may be 
used for JVQ to serve the same purpose. 

Earlier studies on the effects of resonant cavity 



indicat e; 12 ' 13 that both electric and magnetic fields of 
the cavity couple to the Josephson junction since the 
cavity electromagnetic (EM) mode behaves similar to a 
phonon modei^ which interacts with the fluxon. The ef- 
fects of resonant cavity on the fluxon dynamics in LJJ 
stacks^— have been studied both experimentall y 18 ' 19 
and theoretically^ ^— These studies show that when the 
coupling between LJJ and resonant cavity is spatially 
uniform, no force is exerted on the fluxon by the cavity, 
but its dynamics may become modified. These studies 
suggest that the interaction between LJJ and a resonant 
EM wave mode of the cavity promotes^ 3 - collective dy- 
namics of fluxons. The in-phase locking mode of the 
fluxon dynamics is shown to be enhanced^ 3 - by the cavity 
EM mode. 

These studies also suggest that the junction-cavity in- 
teraction may be used to change the qubit property. The 
property of JVQ depends on MQT between two spatially 
separated states of the fluxon. We note that MQT repre- 
sents quantum particle-like collective exciationsi 24 ' 25 As 
semi-classical theories indicate that the MQT rate 2 ^ de- 
pends on the potential barrier height, the JVQ can be 
tuned by adjusting the potential-well for the fluxon. This 
adjustment can be achieved by potential renormalization 
induced by the junction-cavity interaction since this in- 
teraction can strongly affect the fluxon tunneling pro- 
cesses, similar to phonon assisted tunneling in Josephson 
junctions^ We note that a two-level atom interacting 
with a quantized radiation field, described by the Jaynes- 
Cummings model^ is also similar to the JVQ-cavity sys- 
tem that we consider in the present work. The potential 
renormalization for fluxon suggests that the resonant cav- 
ity may be used as a tool for controlling the JVQ prop- 
erty. As the fluxon tunneling processes may be controlled 
externally by tuning either the junction-cavity coupling 
strength or the resonant frequency, the effects of the res- 
onant cavity depend on the nature of the interaction. 
However, the influence of junction-cavity interaction on 
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the MQT rate has not been understood clearly. 

In this paper, we investigate the effects of the junction- 
cavity both on MQT from metastable state and on the 
ground-state energy splitting in a double-well potential. 
We note that, to focus on the interaction between LJJ 
and a single resonant cavity mode, we consider only a 
high-Q c cavity. First, we estimate the MQT rate for 
the fluxon in a single LJJ and for the phase-locked flux- 
ons in a coupled LJJ stack by computing the local and 
non-local contributions. Then, we estimate the effects of 
resonant cavity on the JVQ property by computing the 
ground-state energy splitting. Before proceeding further, 
we outline the main result, (i) The potential barrier for a 
fluxon in the metastable state is not affected by increasing 
neither the junction- cavity interaction nor the resonant 
frequency of the cavity EM mode, (ii) The non-local con- 
tribution to the tunneling rate due to the junction- cavity 
interaction is negligible in the weak-coupling regime, (hi) 
Due to potential renormalization induced by the junction- 
cavity interaction, the potenital barrier height for the 
fluxon trapped in a double-well potential is reduced. This 
reduction leads to increase in the ground-state energy 
splitting for the JVQ with increasing junction- cavity cou- 
pling and resonant frequency. 

The outline of the remainder of the paper is as follows. 
In Sec. II, we describe the LJJ-cavity system by using a 
set of two perturbed sine-Gordon equations. In Sec. Ill, 
the effects of resonant cavity on the fluxon tunneling rate 
from the metastable state in a LJJ are discussed. In Sec. 
IV, we discuss MQT of phase-locked fluxons from the 
metastable state in a vertical stack of two coupled LJJs. 
In Sec. V, the effects of interaction between LJJ and a 
single mode in high-Q c cavity on JVQ are estimated by 
computing the ground-state energy splitting. Finally, we 
summarize the result and conclude in Sec. VI. 



(pi = <fii — 4>i-\). The strength of magnetic induction 
coupling between two LJJs is denoted by S. Here we set 
h = ks = c = 1 for convenience. The perturbation term 
T of for each LJJ which is given by 

= ~dt 9E ~ a 7 t 2 ~ ~ ^ ~ x i) sm V i ( 3 ) 

accounts for the contribution from dissipation (/3), bias 
current (/ = J B / J c ), resonant cavity (<7.e), and microre- 
sistors (e = (J c - J' c )l b / J c Xj). Here x°, J B , J c , J' c , l b 
(<C Xj) and Xj denote the position of microresistors in 
the insulator layer of the i-th junction, the bias current 
density, the critical current density, the modified current 
density, the length of the LJJ in which J c is modified, and 
the Josephson length, respectively. We note that dissi- 
pation, bias currents, resonant cavity and microresistors 
on the phase dynamics lead to different effects. 

We account for the perturbation contribution due to 
resonant cavity by following Tornes and Stroud^ and by 
assuming that the cavity supports a single harmonic os- 
cillator mode which may be represented by the displace- 
ment variable q r as 



d 2 q r co r dq r 2 



(4) 



Here u) r , Q c , and M osc are the dimcnsionless oscillator 
frequency in units of ui p , the cavity quality factor, and 
the "mass" of the oscillator mode, respectively. For sim- 
plicity, we neglect the second term on the left hand side 
of Eq. by assuming that the cavity is non-dissipative 
(i.e., high-Q c cavity). Also, we assume that the cavity 
electric field E is uniform within the junction by consid- 
ering the spatially uniform junction-cavity coupling Qe 
of 



II. COUPLED LONG JOSEPHSON JUNCTIONS 
IN RESONANT CAVITY 

To examine i) one-fluxon tunneling in a single LJJ, 
ii) phase-locked two-fluxon tunneling in a stack of two 
coupled LJJs, and hi) the ground-state energy splitting 
in JVQ, we start with coupled perturbed sine-Gordon 
equations 1 - for describing two LJJs which interact with 
resonant cavity^ 



dt 2 

dt 2 



sin ipi = T\ 
sin ip 2 = T 2 



(1) 

(2) 



where x and t are the dimcnsionless coordinates in units 
of Aj7 _1 (<S) and wr 1 , respectively. Here r y~ 1 (S) = 
\/l — S 2 and uip denotes the plasma frequency. The 
dynamic variable ipi represents the difference between 
the phase 4> of the superconductor order parameter for 
the two superconductor (S) layers i and i — 1 (i.e., 



9E = 'Ye 



M 



4tt 



E • z, 



(5) 



where e<j is the dielectric constant. As we will discuss 
below, the position independent coupling gE does not 
change the fluxon motion directly but yields potential 
renormalization when a microresistor is present. 

To estimate the effects of interaction between LJJ and 
resonant cavity analytically, we consider the weak per- 
turbation T limit. As each perturbation term in Eq. ([3]) 
is small and does not change the form of the kink solution 
within the lowest order approximation^ we describe the 
fluxon motion in terms of the center coordinate q(t). In 
the absence of both the perturbation terms (J 7 = 0) and 
the magnetic induction effect (5 = 0), the fluxon solution 
to Eq. ([T]) is given by 



<Pi(x, t) sa 4 tan 



,n{vi)\x-(n(t)} 



(6) 



in the non-relativistic limit (i.e., « < 1). Here qi{t) = V{t 
denotes the center coordinate for the fluxon, and v is the 
fluxon speed in units of Swihart velocity. Equation ([5]) 
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represents propagation of nonlinear wave as a ballistic 
particle. The perturbation contributions of JF only affect 
the dynamics of fluxon expressed in the q coordinate. 

Wc now describe the fluxon phase dynamics in the cou- 
pled LJJ using the center coordinate qi representation. 
The energy of the fluxon may be seen easily from the 
Euclidean Lagrangian (i.e., r = it), 



C C Q ~t~ C raa g ~l~ Cp er f ~l~ C Q 



Cc 



(7) 



The first three terms for C of Eq. describe the LJJ 
contributions, while the remaining two terms arise from 
the resonant cavity. First, we discuss the LJJ contribu- 
tions to Lagrangian C. The unperturbed part of LJJ is 
described by the Lagrangian C Q given by 



dr 



~&x ) + 2( . :l_a,s > ' 1 



(8) 



The Lagrangian contribution from the magnetic induc- 
tion effect, L mag , is given by 



C mag - S I dx ( ^ j 



(9) 



We note that C ma g accounts for the interaction energy 
Ei n t between two LJJs due to the magnetic induction ef- 
fect. The perturbation contribution to the Lagrangian, 
£pert = £nd + £>d, is expressed as the sum of two terms: 
i) the non-dissipative (C n d) and ii) dissipative (Cd) part. 
The non-dissipative contribution comes from the bias 
currents and microresistors. The non-dissipative La- 
grangian C n d is expressed as the sum of the contributions 
from the bias current (Ctias) and microresistors (C P i n ) 
(i.e., C n d — Cbias + Cpin) ■ The bias current contribution 
Cuas is given by 



(10) 



and the inhomogeneity contribution due to microresistors 
Cpin is given by 



r 



E 



dx eiS(x — x°)(l — cosipi) 



(11) 



We note that C P i n accounts for the fluxon pinning en- 
ergy Epin- These non-dissipative contributions provide 
the bare fluxon potential V(q). On the other hand, the 
dissipative Lagrangian Cd accounts for the interaction 
between the fluxon and environment. The effects of this 
contribution may be described 30 by following Caldeira 
and Leggett and by representing the environment as a 
heat bath. The heat bath is represented as harmonic 
oscillators with generalized momenta Pi and coordinates 
Qi. The dissipation Lagrangian Cd which accounts for 
the coupling between the phase ((f) and oscillator (Qi) 
variables is given by 



C, 



= Jdxj: 



2m 



rriiLuf 



rtutJi 



(12) 



Here, the spectral function Jp(w), 

Jp N = ? E - <*) = (I 3 ) 

Z TfliUJ: 
% L 

is used to reproduce the dissipation effects (fi) in Eq. ((3|) . 
The effects of dissipation on a two-state system has been 
studied extensively by using the spin-boson models In 
the adiabatic approximation, the energy splitting for the 
two-state system is known to be reduced 3 ^ in the dissipa- 
tive environment. However, this result does not^S imply 
that the effects of the interaction between the two-state 
system and a single oscillator, which represents either a 
phonon or quantized radiation field, on the energy split- 
ing is similar. In our discussion below, we neglect the 
dissipation effects by setting (3 = since these effects are 
small at low temperatures, and we focus on the effects 
due to a resonant cavity. 

We now discuss the high-Q c resonant cavity contribu- 
tion to the Lagrangian C of Eq. (J7]) . The resonant cavity 
is modeled by using Lagrangian for a single harmonic os- 
cillator which represents a single EM-mode supported by 
the cavity. The Lagrangian for this single mode oscillator 
C osc is written as 



nsr. Z * 



dq, 
dr 



K 2 

T 9r 



(14) 



where K is the " spring constant" and q r denotes the oscil- 
lator coordinate. We note that the oscillator frequency 
uj r in Eq. Q is given by cu r = (K/Mosc) 1 ^ 2 . The ca- 
pacitive coupling between LJJ and resonant cavity is de- 
scribed by the Lagrangian C coup as 



C 



( dq r 



jdxj: 



dr 



(15) 



Here we assume that the coordinate q r is spatially homo- 
geneous and focus on the effects of the uniform E-field 
in the cavity. We note that the interaction between LJJ 
and resonant cavity yields the non-local effects, similar 
to those from the dissipation term (i.e., (3^0). 

We estimate MQT of fluxon by using the usual semi- 
classical approach 3 ^ of starting with the partition func- 
tion Z for the junction-cavity system 



Z = J T>[Lp]T>[q r } exp{- 



-S[if,q r ]} 



(16) 



where S[ip, q r ] — J drC is the action and C is the La- 
grangian of Eq. (JTJ). By noting that shape distortion 
of the fluxon due to weak perturbation (i.e., small J-) 
is negligible, we may rewrite the partition function Z in 
terms of q(r) and q r {T~) as 



Z = J V[q] J V[q r 



,-S[q,q T ] 



(17) 



Also by noting that the Lagrangian C coup of Eq. (Ti"5)) 
which accounts for the interaction between LJJ and 
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resonant cavity is linear in both coordinates q r and 
tp, we separate the partition function Z into the reso- 
nant cavity and fiuxon contribution by expressing Z — 
Z res Zfi uxon . The resonant cavity (Z res ) and fiuxon 
(-2 fiuxon) contribution to Z are given, respectively, as 

Z res = fT>[q r (u n )] exp{~Sres[q r (^n)]} and Zfl uxon = 

J V[q(T)] exp{— S e ff[q(r)]}. The action for the resonant 
cavity contribution S res [q r ] is given by 



qr,n 



Sres[Qr] — T S 

2Kg E q n u'i 
M osc (u)l + uf) 



M, 



osc II. 2\ 



to 2 ) X 



9r,- 



2irg E q„ 



M osc [wl + rf 



(18) 



where g r ,„ = q,-(u; n ), g n = g(w n ), w n = 2imT is the 
Matsubara frequency, and T is the temperature. The 
action for the fiuxon contribution S e f / [q] is given by 



Seff[q] = / dr 



M e * , 

1=1 

1-S 2 



—2 2 



(Eft)" 



2—1 



1 -S 2 



2 2 

drdr'A-(r-r')V%(r)V ft (r') (19) 



where t& = dqijdr, M e denotes the renormalized fiuxon 
mass 



Mr = M 1 - 



1 2g| 
M l - S 2 



(20) 



due to the spatially uniform junction-cavity interaction 
and the M denotes the rest mass of the fiuxon. The 
mass M e accounts for the renormalization effect of both 
junction-cavity and magnetic induction interaction. The 
bare potential V(q) = V(qi,q2) is given^ by 



z Resonant - 

Uy cavity 




defect 



FIG. 1: A LJJ is shown schematically as an insulator (7) 
layer is sandwiched between two superconductor (S) layers. 
L x and L y denote the dimensions in x— and y— direction, 
respectively. J B denotes the bias current density. The filled 
circle represents microresistor (i.e., pinning center), and the 
dashed box represents resonant cavity. 



After the calcultion , the oscillator coordinate q r in the 
partition function Z of Eq. (fTB|) is decoupled from the 
center coordinate q. This separation allows us to inte- 
grate out the g r -coordinate. Hence, in discussions below, 
we consider the fiuxon contribution Zfi uxon to the par- 
tition function which is described by the action S e ff. 
Using S e ff, we discuss how the junction-cavity interac- 
tion affects both one-fluxon and two-fluxon tunneling in 
LJJs. 



III. MACROSCOPIC QUANTUM TUNNELING 
IN SINGLE JUNCTION 



V(q) 



2c, 



cosh 2 qi 



8S (gi - q 2 ) 
sinh(9i - q 2 ) 



(21) 

Here, the fiuxon potential V(q) includes the effects from 
the three contributions: (i) the potential tilting effect (/), 
(ii) the pinning effect (e) , and (iii) the magnetic induction 
effect (S). The third term in [ ] of Eq. (fL9]) accounts 
for the potential renormalization due to junction-cavity 
interaction. This renormaliztion is similar to that for the 
electronic tunneling process with phonon coupling^ In 
the discussion below, we refer g E = 2ir 2 g E /M osc as the 
strength of junction-cavity interaction. The cavity kernel 
K(t — t') in the third term of Eq. (|T9"]) is given by 



K{r) 



u>$ cosh(w r /2T- 



sinh(w r /2T) 



(22) 



at non-zero temperature T . This term accounts for the 
non-local effect arising from the junction-cavity interac- 
tion. 



We now examine the effects of resonant cavity on MQT 
from the metastable state in a single LJJ obtained by 
implanting a microresistor in the insulator layer and by 
applying the bias current (J B ) as shown in Fig. [TJ The 
dimensions of the junction, compared to the Josephson 
length A j, are chosen so that L x ^> Xj and L y <C Xj. 
These choices are made to enhance the quantum effect 
at low temperatures. We describe MQT of the fiuxon by 
starting with the action S^ff [q] for the LJJ given by 



(It 



e -2 , t/ I \ i -2 2 2 

—q +V s {q)+g E (jJ r q 



g% dTdT'K(r-r r ) q(r)q(T'). (23) 



S effk] 



Here, the action S*jf[q] is obtained from S e ff[q] of Eq. 
(fT9|) . by setting 5 = (i.e., no magnetic induction effect), 
qi = q, and q 2 = 0. Following Caldeira and Leggett, we 
may simplify S^j:[q] by making a usual substitution of 
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V s (x) 




FIG. 2: The fluxon potential V s due to both the bias cur- 
rent density and microresistor in a single LJJ is schematically 
illustrated. 
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FIG. 3: The ratio of the tunneling rates r cQ „(0)/r(0) is plot- 
ted as a function of the junction-cavity coupling strength g% 
to illustrate the size of enhancement. 



q(r)q(r') = [q 2 (r) + q 2 (r')]/2 - [q(r) - q(r'f/2. We 
note that the first two terms of this substitution cancel 
the potential renormalization contribution (i.e., g\u: 2 q 2 
term) arising from the junction-cavity interaction. With 
this cancellation, the action S^ff[q} becomes similar to 
that for the dissipative system^ but the fluxon mass is 
now renormalizcd to 



M e = M 1 



M 



(24) 



and P is replaced by the junction-cavity interaction 
strength (i.e., j3 — > g 2 E ). The renormalized mass M e ac- 
counts for the effects of the uniform E field in the cavity. 
The bare fluxon potential V s (q) is given by 



V s (q)=-2nfq- 



2c 



cosh q 



(25) 



Here the bias current density / = f c — Sf is measured 
in terms of the deviation Sf from the critical value f c — 
4e/(3-\/37r). The potential V a (q) may be approximated 
by a quadratic-cubic potential as shown schematically in 
Fig. [1 The cavity kernel K(t~t') of Eq. (J22D describing 
the non-local effect due to the junction-cavity interaction 
simplifies to 



K{t 



= - Ur \ T -r'\ 



(26) 



in the T = limit. 

The action S^j [q] of Eq. ([23|) indicates that the res- 
onant cavity yields i) fluxon mass renormalization and 
ii) non-local effects. The mass renormalization modifies 
the oscillation frequency about the metastable point, as 
shown in Fig. [5] This change may be easily seen by com- 
puting the oscillation frequency w e at the metastatble 
state (i.e., local minimum) as 



1 d 2 V s (0) 



M, dx 2 



1/2 



l+ 9% 
M 



(27) 



where u is the oscillation frequency at the metastable 
point in the absence of the resonant cavity. The non-local 
contribution due to junction-cavity interaction is similar 
to that for the dissipative system, but to determine the 
size of this contribution more calculation is needed. 

To estimate the size of these two contributions from 
the junction-cavity interaction, we compute the MQT 
rat o 30 i 35 given by 

r co „(0)= A au (0)e- B -(°) (28) 

at T — 0. Here, the prefactor A ca v(0) is given by 

^(0) = V60 We f%^V /2 (29) 



2vr J 

and the bounce exponent B cav (Q) — B 0tCav + SB C 
eludes both the local contribution B Q cav of 



B , c 



(It 



M e . 2 



Vs(q) 



(30) 



and the non-local contribution 5B rnv of 



SB C 



dr / dr'K(r - r') [q(i 

-oo J — OO 



q(r')Y 



(31) 
(0) are 



These two contributions, B 0yCav and 5B cav , to B cav \ 
evaluated explicitly to estimate their size. 

The local contribution B ^ cav may be computed easily 
by approximating V s (q) of Eq. (|2"5)) as a usual quadratic- 
plus-cubic potential of 



V s (x) = V s (q)-V s {q ) 



27V , 2 



(32) 



where x — x/x a , x = q - q , and V — [^-K^S^K^e)] 1 / 2 
is the barrier potential for the fluxon. Here q Q is the 
position of the metastable point and x = 9V3M e ui 2 /32e 
is the escape point as shown in Fig. [2] The evaluation of 
B rav yields 



B„ 



= 2 



dx [2M e V s {x)\ ' = 



1/2 _ 36 V 



5w P 



(33) 



6 



0.010 



SS^y 0.005 



0.000 









0.06 ■ 


/ \ g* = 


0.04 ' 




0.02 


1 / \ . 




"J % V 
1 » 1 




■ / w 




ll /~\ v» 

















0.0 0.2 0.4 0.6 0.8 
CD,. 

FIG. 4: The non-local contribution 8B cav to the bounce ex- 
ponent BcavifS) is plotted as a function of cj r for g% — 0.02 
(solid line), 0.04 (dashed line) and 0.06 (dot-dashed line). 



Using this result, we estimate the local contribution to 
enhancement of the tunneling rate due to the resonant 
cavity. The ratio of the MQT rates, T cav (0)/T(0), is 
given by 



,(0) 



r(o) 



i 



2M 



1+ 72K 

5 co 



(34) 



where T(0) is the tunneling rate in the absence of the res- 
onant cavity (i.e., g% — 0). Equation (f34|) indicates that 
the tunneling rate increases with increasing junction- 
cavity interaction strength g E . In Fig. 02 we plot the 
numerically computed ratio T cav (0)/T(0) as a function 
of g E to illustrate its enhancement in the weak-coupling 
regime (i.e., g 2 E <C 1). The curve indicates that enhance- 
ment of r cat) (0)/r(0) is less than 1%. 

The non-local contribution 5B cav to B cav (0) of Eq. ([25)1 
reduces the tunneling rate T cav (0). The size of this re- 
duction is estimated by evaluating 5B ca v of Eq. (|31j) by 
writing 



SBcav = 



dfdf'e-^ lf - f ' 1 [x(f) - x{f')f 



(35) 

where x(t) — sech (w e r/2). We note that x{t) is the 
solution to the equation of motion for the quadratic-plus- 
cubic potential in the absence of the non-local effect. We 
evaluate Eq. (f35|) and obtain 




JL 



FIG. 5: Two LJJs with a vertical column of two microresistors 
is shown schematically. L x and L y denote the dimensions in 
x— and y— direction, respectively. J B denotes the bias current 
density. The filled circles represent the microresistors. 



To illustrate the cavity frequency dependence, we plot 
SB cav as a function of uj t for g E = 0.02 (solid line), 0.04 
(dashed line), and 0.06 (dot-dashed line) in Fig. £f] The 
curves indicate that SB cav vanishes both in the low and 
high cavity frequency to r limits. Hence, the non-local ef- 
fects on the tunneling rate r ca „(0) is negligible near these 
limits. 



IV. MACROSCOPIC QUANTUM TUNNELING 
IN COUPLED JUNCTIONS 



In this section, we estimate the effects of resonant cav- 
ity on the tunneling rate of the phase- locked fiuxons from 
the metastable state in two coupled LJJs. Here the fiux- 
ons are trapped by the microresistor on each insulator (I) 
layer, shown schematically in Fig. [5] Earlier studies^ in- 
dicate that uncorrelated one-fluxon tunneling is the dom- 
inant process in the absence of resonant cavity. However, 
phase-locking between the fiuxons in two LJJs becomes 
enhanced in the resonant cavity. This enhancement may 
be seen more easily from the effective action S e ff[q] for 
the two coupled LJJs of Eqs. (fH)j) and (|2ip written in 
the rotated coordinates (q+ ,q~) as 



6B C 



_ 2 ( 9V3M e 
29e 1 ~HT~ 



sinh (nuj r /Lu e 



(36) 



The result for SB cav indicates that the non-local contri- 
bution increases almost linearly with g 2 E in the weak- 
coupling regime and has a strong dependence on the 
frequency ui r of the cavity mode. At low cavity fre- 
quencies (u) r <C 1), the non-local contribution varies as 



SB C 



At high cavity frequencies (u> r 3> 1), on 



the other hand, it varies as 5B cav oc w^exp(— 27rw r /o; e ). 



s eff[q] = dr 



M e .2 M -2 rr, s 

-^q\ + —ql + V(q) 



9Wr 2 



J9l_ 
1-S 2 



1-S 2H+ 

dTdr'KiT-T^q+^q+iT 1 ), (37) 



where q± = (qi ± q 2 )/v / 2- The action S e ff[q] indicates 
that the potential for the in-phase mode, (g+, 0), is renor- 
malized by the junction-cavity interaction while the out- 
of-phase mode, (0, q_), is not. Also, the non-local contri- 
bution appears only for the motion in the g+ direction. 
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The bare fluxon potential V(q) = V(q+,q~) of 



V(q +> q-) = -2V2irfq + -- 



8V2Sq- 



-2c 



cosh 



2 q 



sinh \2q^ 



+ 



V2 



cosh' 



\ V2 J 



(38) 



for /i = f'2 = / and c\ — c-i — e, indicates that the 
one-dimensional potential along the (q + , 0) direction (i.e., 
V(q+,0)) corresponding to the in-phase mode becomes 
identical to V s (q) of Eq. (1251 under the transforma- 
tion of 2f — > f, 2e — > e, and q + /^/2 — > q. This sim- 
ilarity reflects that the phase-locked fluxons moving in 
the (q+,0) direction (i.e., qi — (72) behave as a single 
fluxon. However, the one-dimensional potential for the 
out-of-phase mode (i.e., V(0,q_) or along the (0, g_) di- 
rection) behaves as a potential well near the metastable 
point q° = (q+,q°-), determined from the condition 
[dV(q)/dq+] q _ = [dV(q)/dq.} q+ = 0. 

To illustrate these phase-locking modes, we plot the 
potential V(q+,q-) in Fig. © for / = 0.06, e = 0.269, 
and S = —0.05. Here, the value for e and S are chosen 
so that when a vertical stack^l of two interacting JVQs 
are fabricated using coupled LJ Js and microresistors only 
one quantum state is bound on each side of the double- 
well potential. The metastable point q° is denoted by 
the solid circle. The solid lines indicate that the poten- 
tial is metastble for the in-phase mode (i.e., along the 
(q+,0) direction), but it behaves as a well for the out- 
of-phase mode (i.e., along the (0, qJ) direction). These 
curves show that tunneling of the in-phase mode from 
the metastable state is more favorable than that for the 
out-of-phase mode. 

The tunneling rate T cav (Q) from q° can be estimated 
by summing over the contribution from all paths of es- 
cape, but the dominant contribution comes from the 
most probable escape path (MPEP) in which S e ff is 
the minimum.— For the physical parameters chosen in 
Fig. |51 the MPEPs correspond to one-fluxon tunneling, 
indicated by the dashed lines. The MPEPs are deter- 
mined by the two competing energies: (i) the pinning 



energy (£ p 



I Spin I) an d (ii) the magnetic induction 



interaction energy (£ j. 



\E mt \). When £ 

in I > £ v 



the fluxons are not pinned at the microresistor sites but 
maintain a large separation distanced However, when 

£%nt < ^ £', 



pin . 



the one-fluxon tunneling processes are fa- 
vored over the two-fiuxon tunneling processes. 

We now estimate the two-fluxon tunneling rate for the 
in-phase mode. We simplify the calculation by using the 
simialrity between the tunneling of the in-phase mode 
and the one-fluxon tunneling process discussed in Sec. 
III. When the bias current / is less than the critical value 
f c (i.e., f = f c - S f with < S f « f c = 4e/(3\/3V), the 
potential along the path (q + , 0) has the metastable state, 
as illustrated in Fig. [5] The potential V(q+,0) may be 




FIG. 6: The potential VQ(q+,q~.) surface is plotted for e = 
0.269 and 5 = —0.05. The filled circle represents the position 
of the metastable state. The dashed and solid lines denote 
the most probable escape paths (MPEPs) for one-fluxon and 
two-fluxon tunneling, respectively. 



approximated as the quadratic-plus-cubic form of 



271/' 

V(q + ,0)^^(f, 



1+ 



3 ) 



(39) 



where q + = (q + — q + )/q + , q\ is the escape point and 
Vl = 2[d 2 V(q o + ,0)/dql] 3 /3[d 3 V(q + ,0)/dql} 2 denotes the 
potential barrier height for two-fluxon tunneling. We 
note that q + is similar to x Q in Fig. [2j Also, similar to 
the single LJJ, the semiclassically estimated two-fluxon 
tunneling rate of 1^(0) - ^exp[B* OT (0)] at T = 
depends on both the barrier height and oscillation fre- 
quency. The factor J\? cav {fy and bounce exponent Bl av (0) 
are calculated in the same way as in Sec. III. The factor 
-4L„(0) is given by 



Aav(0) « V60 LJ 



BL 



2n 



1/2 



(40) 



The local and non-local contributions to the bounce ex- 
ponents Bl av {0) = B* + SB* av are given by 



R l =2 



dq+^2M e V(q+,0) 



36V? 

5w e 



and 



2f E [9V3M e 



1 



S 2 V 



16e 



sinh (iruj r /oj e ) 



(41) 



(42) 



respectively. The result indicates that the two-fluxon 
tunneling rate r' a „(0) in the cavity is enhanced from that 
r*(0) in its absence. Neglecting the non-local contribu- 
tion, we may write the ratio r' a „(0)/r'(0) as 



rL„(Q) 
r*(o) 



1 



9% 



2M(1 - S 2 



1 + J2Vl 
5 ujo 



(43) 
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This enhancement is similar to the tunneling process 
discussed in Sec. III. The estimated value of r*(0) 
for the Nb-A^O^-Nb-AlaO^-Nb junction is 8.5 x 10 9 
This value is obtained by using the experimental 



valui 



.15.17 



of J c ~2xlO b A/m 2 , X L 



and Up ^90GHz. Also, we chose L y 



90nm, Aj ^25/im, 
^0.2/im to enhance 
the quantum effect and used the experimentally accessi- 
ble valued of e = 0.269, S = -0.05 and 5 f ~5xl0~ 4 . 
On the other hand, the potential V{q +1 q-) along the 
(g + ,0) direction indicates that the two-fluxon tunneling 
rate T* a ^(0) is suppressed from the one-fiuxon tunnel- 
ing rate T° av (Q) along either the q + = q- or q + = — q_ 
direction. This reduction in the tunneling rate is given 
by 



rL„(o) 



r° (o) 

cav V w / 




(44) 



where a D = {[d 2 V(q° + , 0)/dq 2 ]/[d 2 V(q° , Oj/dg 2 ]} 1 / 4 
is a constant of order unity, V° = 
2[d 2 V(q°,0)/dq 2 } 3 /3[d 3 V(q°,0)/dq 3 } 2 is the one-fluxon 
tunneling potential barrier height, V(q, 0) is the fluxon 
potential of Eq. (|21l) along the q + = q_ direction, 
and q° denotes the position of the metastable point 
for one-fluxon tunneling, given by the condition that 
dV(q°, 0)/dq = 0. The ratio r* al) (0)/r° Q1 ,(0) < 1 for the 
potential surface in Fig. [5] reflects that V* 3> V° . 



V. JOSEPHSON VORTEX QUBIT IN 
RESONANT CAVITY 

We now examine the effects of high-Q c resonant cav- 
ity on JVQ. The JVQ may be fabricated by using two 
closely implanted microresistors in the insulator layer of 
the linear LJJ as shown in Fig. [7] As earlier studies^— 
indicate, MQT of fluxon between the spatially separated 
minima of double-well potential leads to splitting of the 
degenerate ground-state energ y 38 ' 39 In this section, we 
estimate the effects of junction-cavity interaction on this 
energy splitting. 

The interaction between the LJJ and resonant cavity 
yields i) fluxon potential renormalization and ii) non- 
local contribution to the action. The effects of these con- 
tributions on the energy splitting may be estimated by 
starting with the action Sfjj for the JVQ given by 



S effil] = 



^q 2 + V Q (q) 



2g% J drdr' if (r 



q(T)q(r') . (45) 



Without loss of generality, we obtain the potential func- 
tion Vq(q) from the double- well potential V(q) of 



V(q) = gW r q 2 



2c 



2c 



cosh [q 



cosh [q 



(46) 



Resonant - 
cavity 




hy///////y///////zm 




defects 



FIG. 7: A LJJ with two microresistors, representing a Joseph- 
son vortex qubit, in a resonant cavity is shown schematically. 
The separation distance between the microresistors is denoted 
by t. The filled circles and dashed box represent the microre- 
sistors and resonant cavity, respectively. 



8,= 8 ? =S 



g;^0 / g E"° 




FIG. 8: A schematic diagram of a double- well potential Vg(q) 
due to the two microresistors in the insulator layer of the LJJ 
is shown to illustrate the renormalization of Vg(q). The solid 
and dashed lines represent the potential Vg(q) in the absence 
and in the presence of the resonant cavity, respectively. 



where £ denotes the separation distance between the two 
microresistors. Here, we have added a constant energy 
E Q term to V(q) (i.e., V Q (q) = V(q)+E Q ) so that V Q (q) 
vanishes at the potenital minima. Here, the potential 
Vq((?) may be characterized by the position of the two 
minima and the potential barrier height. In the discus- 
sion below, we do not make the usual substitution of 
q{T)q(r>) = [q 2 (r) + q 2 {r')]/2 - [g(r) - q(r')} 2 /2 used in 
Sec. III. This approach allows us to elucidate the origin 
of the changes in the energy splitting due to the junction- 
cavity interaction. 

In the absence of the resonant cavity (i.e., g% — 0), 
the double- well structure for Vq(<j) with the separation 
distance £ > £ « 1.317 is shown schematically in Fig. [5] 
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FIG. 9: The ratio of the potential barrier height V ,cav/V is 
plotted as a function of the junction-cavity coupling strength 
g% for Lo r = 0.50 (dot-dashed line), 0.70 (dashed line), and 
0.90 (solid line) to illustrate the suppression in the cavity. 



FIG. 10: The shift 5 in the position of the potential minima is 
plotted as a function of the junction-cavity coupling strength 
g% for L0 r = 0.50 (dot-dashed line), 0.70 (dashed line), and 
0.90 (solid line). 



as the solid line. The two potential minima are located 
at q = ±q /2 where q Q is determined from 



cosh q Q 



sinh 2 i - 1 
cosh^ 



(47) 



The energy shift Eq, representing a constant of motion, 
is given by 



E Q = -2e- 



cosh t 
cosh 2 ^- 1 



(48) 



Also, the potential barrier height V between the two 
minima (i.e., q = ±q a /2) is given by 



V = 2e 



/cosbJ- 2 
V sinh^ 



(49) 



We note that these quantities change in the resonant cav- 
ity, as shown schematically by the dashed line in Fig. [H 
In the resonant cavity (i.e., g 2 E ^ 0), on the other hand, 
the JVQ potential Vq(q) acquires an additional g E uj 2 q 2 
term in Eq. (j46]). This term arises from the coupling 
between the oscillator coordinate q r and the center co- 
ordinate q in the coupling Lagrangian C coup of Eq. (fT5)) 
and accounts for potential renormalization. The main 
renormalization effects are the following: i) the barrier 
potential height is reduced, ii) the position of the po- 
tential minima become closer together, and iii) the os- 
cillation frequency at the potential minima is modified. 
These effects become amplified with increasing junction- 
cavity interaction strength (g%) and resonant frequency 
(oj r ). 

The effects of the junction-cavity interaction on the 
potential barrier height V 0jCav may be estimated straight- 
forwardly. In Fig. IH1 we plot the numerically computed 
ratio V 0tCav /V as a function of g 2 E to illustrate the de- 
pendence on the junction-cavity interaction. The curves 



for uj r — 0.50 (dot-dashed line), 0.70 (dashed line) and 
0.90 (solid line) indicate that the barrier potential height 
decreases with increasing g E and uj r . Also, the curves 
indicate that the ratio decreases linearly in the weak cou- 
pling regime. To leading order in g 2 El the potential barrier 
height V 0<cav estimated from the renormalizcd potential 
V{q) of Eq. (H6T) is given by 

(50) 



-2 2 9 



This decrease in the potential barrier height leads to the 
increase in the ground-state energy splitting. 

Another important effect of the resonant cavity is the 
shift S a in the position of potential minima. As the poten- 
tial barrier height is reduced, the position of the potential 
minima are closer together. The shift S a from the initial 
position of q = ±q /2 is given by 



5n = 



2 cosh 2 nanh 2 ^ 



cosh2£- 7 



(51) 



Here, we obtained 5 by imposing the condition 
[dV{q)/dq\ qHqo/2)± = 0, where (q /2) ± = ±[{q /2) - 5 ] 
denotes the new potential minima. This shift S modifies 
the constant of motion Eq. The new value for Eq may be 
obtained from the condition [dq(r) / dr]tg /2) ± — 0, noting 
that the fluxon is initially located at the bottom of either 
side of the double- well potential so that Vg((g /2)±) = 0. 
We plot the numerically computed shift 5 a as a func- 
tion of g 2 E in Fig. [TU] for uj r — 0.50 (dot-dashed line), 
0.70 (dashed line), and 0.90 (solid line) to illustrate the 
amount of this shift in the weak-coupling regime. The 
curves indicate that 5 increases with g 2 E and with ui r , 
reflecting potential renormalization. 

The resonent cavity also modifies the oscillation fre- 
quency w e at the potential minima. The modified fre- 
quency u e is given by 



1 + 

M 



1 



(52) 
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where lo is the frequency in the absence of resonant cav- 
ity and T = 6q a sinh 2q tanh £/e(cosh 2 £ — 4) sinh 2 1. 

We now combine these effects together and estimate 
the ground-state energy splitting^ A cav by using the 
action Sfj^[q] of Eq. (j4"5"j) and by using the standard 

method of summing over the "instanton" trajectories,— 
By following Weiss and coworkers^ we compute the one- 
bounce contribution to the partition function Zfi uxon , 
assuming that the fluxon is initially pinned at one of the 
potential minima. We write the partition function as 



-'fluxon 



(53) 



i=0 



where Zi denotes the i-bounce contribution. Here the 
bounce is an instanton-anti-instanton pair. To estimate 
A cm , we compute both the saddle-point {Z ) and the 
one-bounce (Z\) contribution to Zfi uxon by noting that 
Z\ may be expressed as 



„ Zq ( A cav t 
Zl - — 



(54) 



where = 1/T. For the contribution Zq, we assume that 
the fluxon is initially confined at q = (g Q /2)_ and obtain 



(55) 




where the eigenvalues A° are determined from 



q 



Me dt + V Q (-^ + S) q° n {r) 
0/2 



+Airg% / K(t - r') q° n (r') = «(r).(56) 

J —9/2 

Here d 2 = d 2 /dr 2 , Vg(q) = d 2 Vq(q) / 'dq 2 , and the cavity 
kernel K(t — t') = (w^/2) exp[— w r |r — r'|] accounts for 
the non-local effect. 

For the one-bounce contribution Z\ to Zfi uxon , we sep- 
arate the center coordinate q(r) into two parts as 



q(r) = q{r) + £ CnQnir) , 



n=0 



where q(r) describes a bounce- like trajectory and the re- 
maining terms describe the arbitrary paths about this 
bounce-like trajectory. This separation of q(r) may be 
used to write the action S%j [q] as 

OO -j 

Sf ff [q(r)} = S c B ^(q(r)) + £ -\ n c 2 n . (58) 

71=0 

Here S c q\ accounts for the one-bounce-like trajectory in 
the resonant cavity. We choose q n {j) of Eq. (|5"7| so that 



the eigenfunctions of the second variational derivative of 
S®jj:[q\ at q and the eigenvalues A„ are determined from 



-M e d 2 + V Q (q) q n (r) 



^9 2 E 



0/2 



K(t-t') q n {T') = A n g„(r). (59) 



-9/2 



We note that the first two eigenvalues, Ao and Ai, need 
to be separated from the rest because Ao < and Ai = 
while the other eigenvalues are positive. The one-bounce 
contribution (Zi) may be expressed as 



Zi = N U 



dc„ 



2tt 



(60) 



where N is a normalization constant. With the separa- 
tion of the first two eigenvalues (i.e., Ao < and Ai = 0) 
from the others, we write the one-bounce contribution to 
the partition function as 



Zx 



Z 9 
2tt 



■ ,e 
Jo 



n, l= 2 A ™ 



-9/2 V« T 1 



1/2 r 



-9/2 \ dT 



1/2 



-,1/2 



(61) 



We now need to evaluate Z\ of Eq. (|6Tj) to estimate 
A cav . Using Eq. (|54]l . we write the ground-state energy 
splitting A cav as 



A c(liV 



2u e / 



RcavL 



■cav J - / cav^ 



1/2 



where the dimensionless factors R ca v and L cav are 

1/2 



Re 



1 



noo , 

Yln=2 *n 



(62) 



(63) 



and 



dr 



dq 
dri 



1/2 



dr' 



dq 
d? 



1/2 



(57) respectively. The exponent Sjf \ is given by 



Qcav 
D B,1 



0/2 



dr 



-9/2 



Me f dqjr) 
2 V dr 



VqGz) 



(64) 



(65) 



This exponent accounts for the contribution from the two 
transversal of the potential barrier. We note that the ex- 
ponent Sg'\ of Eq. (|65|) does not contain the non-local 
contribution, as in Eq. (|2"8")l . because this contribution is 
already included in the calculation of Z\ (see Eq. ([60)) '). 
We now compute R cav , L cav and S'g\ , separately, to de- 
termine the ground-state energy splitting A cav . To focus 
on the effects due to the junction-cavity interaction, we 
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8e 

FIG. 11: The numerically computed ratio of the dimensionless 
factor R C av/R is plotted as a function of the junction-cavity 
coupling strength g% for ui r = 0.60 (dot-dashed line), 0.75 
(dashed line), and 0.90 (solid line). 




FIG. 12: The numerically computed ratio of L cav /L is plotted 
as a function of the junction-cavity coupling strength g% for 
uj r = 0.60 (dot-dashed line), 0.75 (dashed line), and 0.90 (solid 
line) to illustrate the enhancement. 



present the details of the calculation for R cav and L cav 
in Appendix A and B, respectively, and discuss the de- 
pendence of these factors on the junction-cavity coupling 
strength g\. 

The dimensionless factor R cav in the weak-coupling 
regime is given by 



R r 



-2 2 

ng E ^ 



Xr 



2Muj 2 (a 



(66) 



where Xr — + ^5uj 2 uj + \2uj r uj 2 — 2u; 3 . Equation 
(|66|) yields the value R ca v = 2 in the absence of resonant 
cavity (i.e., g\ — 0). 39 In Fig. [TTJ we plot the numerically 
computed ratio R cav /R as a function of g 2 E for uj r = 0.60 
(dot-dashed line), 0.75 (dashed line), and 0.90 (solid line) 
to illustrate enhancement of R due to resonant cavity. 
The curves indicate that R cav /R increases from 1 almost 
linearly with increasing g E and uj r . 

For the dimensionless factor L cav , we evaluate the in- 
tegral of Eq. by expanding the function Q(t) which 
accounts for the non-local contribution to the bounce-like 
trajectory as a power series. (See Appendix B.) In the 
weak-coupling regime (i.e., g 2 E <C 1), we obtain 



smaller value than that for R cav - To illustrate this de- 
viation, we numerically compute L cav and plot the ratio 
Lcav/L in Fig. [T2]as a function of g% for u> r = 0.60 (dot- 
dashed line), 0.75 (dashed line), and 0.90 (solid line). The 
curves show nonlinear enhancement of the dimensionless 
factor L cav for much smaller value of g E than that for 
R cav shown in Fig. 111! 

Finally, we estimate the effects of junction-cavity in- 
teraction on Sg™. The action of Eq. §5$ for the 
bounce-like trajectory is given by 



(9o/2)h 



(<3o/2)_ 



dq^2M e V Q {q) . (68) 



The integral of Eq. (|68| is evaluated in the same way as 
that for Lcav (see Appendix B). Again, we simplify the 
calculation by writing Vq as a power series in q and then 
expand \/Vq(<?) in powers of g E as 



V M 



V2M e q 
~1 



1 - 



2eq 2 



V 



h ^ Ah2 2 
h + —q 



-9% 



2 2 

oj a 



8 2V 



^ U X \ 2 

6i H — b 2 q 

2V \ 3 y 



, (69) 



V M [A +g 2 E (B Q + B 2 q 2 +B 4 q 4 )} 



(67) 



by retaining the leading order contribution (in g%). Here 
V M = qoVZMVo, A = l-q 2 {2eb 1 /W )-q*(4eb 2 /lW ), 
B = -(e + 8b 3 uj 2 )/8e, B 2 = [he + 2(66i6 3 - 
1)oj 2 + 2 7 rdio; 3 ]/6V , and B 4 = [b 2 e + 20b 2 b 3 u} 2 + e) + 
3TTd3L0^]/15V o . The frequency independent constants bi 
are b x = (cosh£- 2)sech 4 (^/2), b 2 = (cosh 21- 26 cosh 1 + 
33)/(cosh£+l) 3 , and b 3 = (sinh^tanh£) 2 /(cosh2^ - 7). 
Equation ([67]) indicates that L cav in the resonant cavity 
is larger than L = VmA in its absence. However, due 
to the functional form of L ca v, the enhancement of L cav 
from L deviates from the linear dependence on g E at a 



where b\ = b± + (lu 2 /e). Using this series expansion for 
\/Vq, we evaluate Eq. (j68| and obtain S B a \ to the lead- 
ing order in g E as 



Qcav 



AVm[A + g%(B + B 2 q 2 + B iq A )] 



(70) 



where B 2 = B 2 — (ndiUJ^/3V ) and B4 = B4 — 
(3ird 3 ujf/15V ). The action \ in the presence of cavity 
is reduced from that in its absence (i.e., S B a \ < Sb,i)- 
To illustrate this suppression of the ratio, we plot the 
numerically computed ratio S B a \/SB,i as a function of 
g 2 E for oj r — 0.60 (dot-dashed line), 0.75 (dashed line), 
and 0.90 (solid line) in Fig. [T3] The curves indicate that 
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FIG. 13: The numerically computed ratio of the action 
Sj3i/Sb,i is plotted as a function of g% for oj r = 0.60 (dot- 
dashed line), 0.75 (dashed line), and 0.90 (solid line) to illus- 
trate that the one-bounce-like action is reduced. 



FIG. 14: The numerically computed ratio of A ca «/A is plot- 
ted as a function of the junction-cavity coupling strength g% 
for uj r = 0.60 (dot-dashed line), 0.75 (dashed line) and 0.90 
(solid line) to illustrate the enhancement in resonant cavity. 



in the one-bounce contribution to the action decreases 
almost linearly with g E in the weak-coupling region as 
indicated by Eq. (l70l) . This reduction reflects that the 
potential barrier height is reduced (see Fig. [5]) and the 
potential minima become closer together (see Fig. [T0|) 
with increasing junction-cavity interaction strength. 

We now combine the effects of resonant cavity on R cav , 
L cav and Sb,i together and estimate the enhancement of 
the gound-state energy splitting A cav from A. Here, A 
denotes the energy splitting in the absence of resonant 
cavity given by 



■a \ 1/2 

,4, |) e- 



(71) 



where A = UI^Lo K./ Il^Li ^«] 1/2 and S a denotes the 
action integral. In the weak-coupling regime, the ratio 
A ca „/A to the leading order in g E is given by 



^■cav _ 1 g E L u[ 

A ~ M\ uj 2 
M 



ttX 



It 



+ 



8(u; r +w ) 3 _ 
(Bo + B 2 q 2 Q + B 4 q 4 a ) (72) 



2A 

2MV M (B + B 2 q 2 o + B 4 q 4 o ) 



The result indicates that A cav is enhanced with increas- 
ing g E and w r . To illustrate this enhancement, we nu- 
merically compute and plot A ca „/A as a function of g E 
for uj r — 0.60 (dot-dashed line), 0.75 (dashed line) and 
0.90 (solid line) in Fig. The curves show that A cav /A 
increases roughly linearly with g E from g E = to 0.02. 
However, the deviation from this linear behavior becomes 
noticeable for g E > 0.02. Also, A cav /A increases signifi- 
cantly from 1 at g E — in the weak-coupling regime. We 
note that the corresponding changes in the ratio R cav /R, 
Lcav/L and S c e v 1 /Sb,i over the same range of g E are 



less significant. For instance, A cav /A for uj r = 0.90 in- 
creases from 1.0 to 1.45 for the increase of g\ from 0.0 
to 0.015. Over the same range of g%, Rcav/R, L cav /L 
and S c B a ^/S B ,i change from 1.0 to 1.05, from 1.0 to 1.29, 
and from 1.0 to 0.96, respectively. The notable increase 
in A cau /A compared to R cav /R, L cav /L and S^/S B ,i 
reflects that A is small. 4 Hence, A cav depends sensitively 
on the variation of the exponent S B a \. 



VI. SUMMARY AND CONCLUSION 

In summary, we investigated the effects of high-Q c res- 
onant cavity on MQT of fluxon from metastable state in 
a single LJJ and in a stack of two coupled LJJs. Also, we 
estimated the ground-state energy splitting for fluxon in 
a double-well potential. We find that both the tunneling 
rate and the ground-state energy splitting are increased 
in the resonant cavity. However, the amount of these in- 
creases is significantly different. For MQT of the fluxon, 
the tunneling rate increases due to the renormalization of 
fluxon mass, but negligible in the weak-coupling regime. 
On the other hand, the increase in the ground-state en- 
ergy splitting is due to potential renormalization, but 
this increase can become significant with increasing g E 
as shown in Fig. 1141 This energy splitting enhancement 
is consistent with the result of increase in the energy sep- 
aration due to the interaction between a two-level system 
and a quantized radiation field, described by the Jaynes- 
Cummings (JC) models Moreover, the consistency- 4 ^ 
between the result of the present work and that of the 
JC model indicates that the effective Hamiltonian for the 
JVQ-cavity system may be similar to the JC model. 

The effects due to i) interaction between the JVQ and 
a dissipative environment and ii) the losses resulting from 
a low-Q cavity are neglected in the present work. These 
dissipative effects are expected to be present in real sys- 
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terns and may be accounted by using an effective spectral 
density which characterizes the form of dissipation^ In- 
clusion of both the dissipative environment and cavity 
losses may reduce the size of increase in the ground-state 
energy splitting and may lead to decrease in the energy 
spliting when the dissipative effects become strong, as in- 
dicated by the analysis of dissipative two-state systems^ 
However, these dissipation contributions do not reverse 
the effects due to the potential renormalization com- 
pletely in weakly dissipative systems. 

Enhancement of ground-state energy splitting due to 
the junction-cavity interaction may have an important 
consequence for the decoherence time of JVQ in the res- 
onant cavity. Earlier study^ of the JVQ decoherence 
time by Kim, Dhungana and Park indicates that the in- 
crease in the decoherence time in noisy environment (i.e., 
T£ mse ) is correlated with the increasing ground-state en- 
ergy splitting A. This suggests that, as A may be tuned 
by adjusting the strength of junction-cavity interaction, 
the resonant cavity may be used to control the prop- 
erty of JVQ. For instance, the decoherence time T£ mse 
may be increased by increasing the strength of interac- 
tion between fluxon and cavity EM mode. Also, due to 
the similarities between a cavity EM mode and an optical 
phonon mode, the interaction between fluxon and optical 
phonons in the LJJ may affect the decoherence time. 

Another important property of JVQs is entanglement 
between the qubits. As our result suggests that the deco- 
herence time for JVQ can be increased by increasing the 
strength of junction-cavity interaction, the resonant cav- 
ity may also be useful for tuning the level of entanglement 
between the JVQs. Our study suggests that the present 
approach for JVQs is similar to the microwave cavity ap- 
proach used for the other superconductor qubits4^ The 
effective Hamiltonian for the multiple JVQs in a resonant 
cavity may resemble the Tavis-Cummings model^ which 
is the extension of the JC model to the case of multiple 
qubits. This similarity may be exploited by using the 
resonant cavity to control the level of concurrence^ for 
JVQs since the junction-cavity interaction may also pro- 
mote entanglement. Hence, the effects of resonant cavity 
on entanglement between the interacting JVQs would be 
an interesting area for further study. 



where 8± (A) denotes the phase shift due to the scattering 
potential U. This phase shift may be expressed as 



<5±(A) = cot 



U- 1 - g' x (0) T g' x (r s ) 
g'm±g'{{r s ) 



(74) 



where t s = —8/2, and g' x (r) and g x {r) denote the 
real and imaginary part of the Green's function (i.e., 
<7a(t) = g'\{r) + ig'xir)). The phase shift <5±(A) due to 
the scattering from the net potential difference of 



VMq)-VMq /2)- 



U 



8 



consists of two Dirac 8- functions at t = ±r s /2. 
strength of the scattering potential U is given by 

U- 1 = <?o(0) - <?o(t s ) , 



(75) 
The 

(76) 



where go (r) is the Green's function for the eigenvalue 
A = 0. The Green's function g x (r) is written as 



5a (t) 



d'jj 



_ x 2tt M e [uj 2 + CM + wf] - A - iS 



(77) 



Here the effects of the resonant cavity are accounted for 
via M e , uj e and The function C(w), obtained from 

the cavity kernel K{t) of Eq. (f2"fij). 



CM 



Aitg 



Me. 



(78) 



reflects that the resonant cavity supports a single-mode 
with frequency u) r . Using the function C( w )i we write the 
real part of the Green's function as g' x (t) — g' x , (r) + 
g'x -( r )i where 



-1 



AM e u>\ 



1 ± 



Ha 



sinu;A,±T , (79) 



1/2 



= [(A/M e ) - ^ ~ ^ 2 ]/2, 



and ojI x = {[(A/M e ) - ^ + oj 2 } 2 - (Wwgl/M^} 1 / 2 . 
On the other hand, we write the imaginary part of the 
Green's function as <7 a (t) = g' x ' , (t) + g x _(t), where 



The authors would like to thank W. Schwalm and K.- 
S. Park for helpful discussions and I. D. O'Bryant for 
assisting with part of the numerical calculation. 



APPENDIX A: 

CALCULATION OF R cav 

For convenience, the dimensionless factor R cav of Eq. 
(|63|) is estimated in the continuum limit. In this limit, 
we may write R cav as 



Rcav — exp 



dX 



[S + (X) + 8-(X)} 



(73) 



9'L ± {r) 



1 



4M e w A . : 



1± 



n 2 



- u 



1,A 



2^,A 



COSWA,±T . (80) 



We note that the phase shift 5±(A) has both slowly vary- 
ing and rapidly oscillating contributions. For an ex- 
tended bounce (i.e., cu e T s ^> 1), the rapidly oscillating 
terms become negligible compared to the non-oscillating 
terms. 

The factor R cav of Eq. (|73|) may be simplified by us- 
ing the substitution A = M e uj 2 (l + p 2 ), where p is a 
dimensionless momentum variable. With this change of 
variable, we write R cav as 



R cav = exp i - 

7T 



p dp 
l+p 2 



[S+(p) + S-(p)] 



(81) 
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The factor R cav of Eq. (1811) may be further simplified 
by neglecting the rapidly oscillating contributions in the 
phase shift <5±(A) of Eq. ([74]) . Neglecting these oscilla- 
tory contributions, we approximate S±(p) to a simpler 
form 5{p) and write the factor R cav as 



R c 



exp 



p dp 



Sip) 



1+p 

The simplified phase shift S(p) is given by 



(82) 



S(p) = cot 



where the scattering potential strength U is given by 




(84) 



and Wo — {ui 2 + uif )/lu% . We note that ivifl and uj2,q 
are obtained from u>i t \ and u>2,\ of Eq. ()79|) for the eigen- 
value A = 0, respectively. The real and imaginary part 
of the Green's function are given, respectively, by 



and 



i 



2Mpui 



8Mp 2 u 2 {p 2 ujl+u) 2 f 



(90) 



where X g = 4oj 2 [ttoj 2 (3p 2 ojI + oj 2 ) - 2Tp 2 (p 2 uj 2 + u 2 ) 2 ] + 
Mp 2 uj 2 (p 2 uj 2 + lu 2 ) 2 . Now, we use Eqs. §55} - and 
rewrite the simplified phase shift 5(p) of Eq. ([53"} as 



5 (p) = cot 1 p 



2Mp(l+p 2 )X u _ 



(91) 



(83) where X p = -(w*w? + 2w e ^ + u*) + p 2 {2uj 5 + Ljfcj r 



4u}^uj 2 - 3w>j? - 8u> u)* - 4cj- 5 ) - P 4 {16oj^uj 2 + 8uj 2 oj?) - 
p 6 {8to 5 + Aw A w r ) and X^ = uj 2 (uj + uj r ) 2 (p 2 u 2 +co 2 ) 2 . 
Finally, we substitute 8(p) of Eq. (|9ip into R cav of Eq. 
(|82[) and evaluate the integral to obtain 



Rc 



ng 2 E uj 2 



Xr 



2Mlu 2 (u> r + 0J o ) 3 ' 



(92) 



where Xr = 5ujf + 15uj 2 u} + 12uJrUJg-2u}g. Equation ([92} 
yields R ca v — 2 in the absence of the resonant cavity (i.e., 
g% = 0) as expected^ 



and 



l 



4M R 



(85) 



4M R 



(86) 



where W p = (u) 2 +lj 2 p )/uj 2 „. We note that u>i tP and U2, P 
are obtained from and LU2.\ of Eq. (|79[) . respectively, 
by setting A = M e uif(l +p 2 ). 

We now compute R ca v to the leading order in g E to 
account for the effects of resonant cavity in the weak 
coupling regime (i.e., g 2 E <C 1). For this calculation, we 
write the renormalized mass of the fluxon as M e = M — 
2g E and express the oscillation frequency u e as 



^ = uj 2 



1 



Ml 
M 



1 



(87) 



APPENDIX B: 

CALCULATION OF L cav 

The factor L cav of Eq. (|64p may be estimated by deter- 
mining the bounce- like trajectories <?(t). The trajectories 
obey the equation of motion given by 



M, 



d 2 q(r) dV Q {q) 



dr 2 



dq 
^f E 



dr'K(T - r')g(r) = .(93) 



We rewrite the equation of motion in a convenient form 
by integrating Eq. (|93|) by parts and obtain 



Me (dq 
2 [dr 



Vq(q) 



+4irg 



dr'K(T - T r )q(T)q(T r ) = 0. (94) 



Also we rewrite the strength of the potential U as 
1 ^9 2 E ^r ( - x 



2Mu) M 2 (u 2 -lu 2 ) 2 



1 



327rwi?a;3 



(88) 



where X u = Mw 2 [u 2 - uj 2 ) 2 - 8[cj 2 T(u; 2 - to 2 r ) 2 + 
2irujf(uj 2 — 3uj 2 )]. By combining these expressions to- 
gether, we rewrite the real and imaginary part of the 
Green's function of Eqs. (1551 and ([55} , respectively, as 



Using this result, we write the factor L cav as 

2 



L, 



2 J \dr 



K<Zo/2) + 
'(9o/2)- 

where q = q(r) and 



dqsJv Q (q) + 27: g 2 E u 2 r qQ(T) (95) 
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q(t) 





1 


Y 


J 


-2+5 



Q(x) 



FIG. 15: Similarity between the function Q(r) of Eq. 
and the instanton solution q(r) representing the trajectory 
of the fluxon from one potential minimum to the other via 
tunneling is illustrated schematically. 



Here, the non-local contribution due to resonant cavity 
is accounted for by Q(t). As discussed in Appendix C, 
the function Q(t) is similar to q(r). By exploiting this 
similarity, we expand Q(t) in a power series as 



Q(r) 



n=0 



hn+ll 



2)H 



(97) 



where <i 2 n+i is the expansion coefficients (see Appendix 
C). The power series expansion for Q(t) allows us to 
evaluate the factor L cav straightforwardly. By using this 
power series expansion, we evaluate the integral of Eq. 
(|9"3")) in the weak-coupling regime (i.e., g E <C 1) and ob- 
tain the factor L cav to the leading order in g E as 



L C av ~ Vm[A + Qe (Bo 



(98) 



where V M = q ^2MV , A Q = 1 - ^(2e&i/3V ) - 
q 4 (4eb 2 /15V ), B Q = -(e + 86 3 ^)/8e, B 2 = [he + 
2(6b 1 b 3 -l)cu? + 2Trd 1 U}3]/6V , and B 4 = [b 2 e + 20b 2 b 3 uj 2 + 
e) + 3ird 3 Lj3yi5V . The frequency independent con- 
stants bi are given by b\ — (cosh^ — 2)sech 4 (£/2), 
6 2 = (cosh2^ - 26 cosh ^ + 33)/(cosh£ + l) 3 , and b 3 = 
(sinh I tanh I) 2 / (cosh 2£-7). 



APPENDIX C: 

POWER SERIES EXPANSION OF Q(r) 

The numerically computed function Q(t) of Eq. (|9"o) 
indicates that Q(t) is similar to the functional form of the 
bounce-like trajectory q(r). This similarity suggests that 
Q(t) is a scaled function of q(r) as shown schematically 
in Fig. 1151 In this case, we may express the function 



Q{t) as a power series in q(r) as 

oo 

Q(r) = Y,d2n+iq 2n+1 (r) , 



(99) 



n=0 



where d 2n +\ denotes the coefficient for this power series 
expansion. We compute the coefficients d 2n+ i by starting 
with a series expansion of q(r) in r as 



j(T)=Y j a2n+lT 2n+1 , 



(100) 



n=0 



noting that the instanton solution q(j) is an odd function 
of t. Here, the coefficient d 2n +i is obtained by following 
the five steps as discussed below. First, we write the 
bounce-like trajectory q in the absence of resonant cavity. 
This trajectory q may be expressed as 



q 



-b\T + b 2 tanh 1 (63 tanh g), 



(101) 



where the constants b\ = 2y/e/M coth^, b 2 — (cosh2g + 
cosh^)/ sinh2g D , and 63 = cothg depend on the param- 
eters I and e. Second, we expand the right hand side of 
Eq. (llOip as a power series in q as 



-b 1 T+b 2 b 3 q 1 



2-5bj + 3b* 4 
15 q 



(102) 

Here, we find the coefficients a 2n+i by substituting the 
series expansion for q(r) of Eq. (I100|) into Eq. (|102l) . 
The first three coefficients are given by 



Ol 



«3 



«5 



1>1 



b 2 b 3 - 1 
blb 2 b 3 (l - b 2 3 ) 
3(6263 -I) 4 ' 

6f6 2 6 3 (l - 6§)[6 2 6 3 (3 - 26|) + (2 - 36j)] 
15(6 2 6 3 - l) 7 



Third, we use Eqs. (J26]) and (ITOO)) to evaluate Q(r) of 
Eq. (|TJ51) explicitly as 



00 />oo 

Q{t) = ^ a 2n+1 / dr' 



r\r— t'| I2n+X 



(103) 



Fourth, we evaluate the integrals of Eq. (|103p and write 
Q(t) in a power series in r as 



Q(r) » - 



r ai 



6« 3 

,.,2 



120os 



t 3 I a 3 



20a 5 \ 



+ r°a 5 



(104) 



Finally, we use the power series expansion for q(r) of Eq. 
([TOOl) and rewrite Q(r) of Eq. flM]) as 



Q(r) = r(diai) + r 3 (d ia3 + d 3 a\) 

+ r 5 (d 1 a 5 + 2>d 3 a\a 3 + d 5 a\) 



(105) 
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This series expansion allows us to obtain the expansion 
coefficients d 2n +i by comparing the power series Q(t) of 
Eqs. (|104|) and (|105j) . The first three expansion coeffi- 
cients, tfen+l) are the following: 



dn = 



12 



7a 7 



lla 3 a 5 3a| 



di 
d 3 



2 
4 



_6_03 



( Was 
V of 



3a§ 



60 



7a7 



a 3 a 5 



In Sec. V, we use these expansion coefficients to esti- 
mate the dimensionless factor L cav and the one-bounce 
contribution to the action (i.e., S^T). 
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